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We present quantum-theoretical studies of collisions between an open-shell S-state atom and a 2 II-state molecule in the presence 

r*| of a magnetic field. We analyze the collisional Hamiltonian and discuss possible mechanisms for inelastic collisions in such sys- 

Qh terns. The theory is applied to the collisions of the nitrogen atom ( 4 S) with the OH molecule, with both collision partners initially 

"i in fully spin-stretched (magnetically trappable) states, assuming that the interaction takes place exclusively on the two high-spin 

G (quintet) potential energy surfaces. The surfaces for the quintet states are obtained from spin-unrestricted coupled-cluster cal- 

5 culations with single, double, and noniterative triple excitations. We find substantial inelasticity, arising from strong couplings 

*T\ due to the anisotropy of the interaction potential and the anisotropic spin-spin dipolar interaction. The mechanism involving the 

• dipolar interaction dominates for small magnetic field strengths and ultralow collision energies, while the mechanism involving 

Q the potential anisotropy prevails when the field strength is larger (above 100 G) or the collision energy is higher (above 1 mK). 

The numerical results suggest that sympathetic cooling of magnetically trapped OH by collisions with ultracold N atoms will not 

>~>be successful at higher temperatures. 
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The first experimental realization of Bose-Einstein conden- 
sation in a dilute gas in 1995 ' opened up a novel and fast- 
growing field of research on cold and ultracold matter. At 
temperatures below about 10 -6 K, novel properties emerge 
in which the quantum nature of atoms and molecules is cru- 
cial. Although the original experiments involved quantum- 
degenerate states in atomic systems, it was soon realised that 
molecules, especially those with a permanent dipole moment, 
offer an additional range of applications in physics and chem- 
istry. These include development of new frequency standards, 
tests of fundamental physical concepts such as parity and 
time -reversal violation 2,3 , spectroscopic measurements of un- 
precedented accuracy 4,5 , quantum information processing, 6,7 , 
and control of chemical reactions with state-selected reagents 
and products^—. 

In contrast to atoms, which nowadays can be cooled rela- 
tively easily by laser Doppler cooling and evaporative cool- 
ing—, molecules are incomparably more challenging because 
of their complicated internal structure. Two main classes of 
methods have been established to produce cold molecules: di- 
rect methods, in which molecules are cooled from high tem- 
perature by means of a buffer gas or external fields, and indi- 
rect methods, in which cold molecules are formed from pre- 
cooled atoms by photoassociation or magnetoassociation. 
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Indirect methods can now produce ground-state molecules 
at temperatures below 1 jSKM ~— . It has been shown recently 
that, for KRb, the rates of chemical reactions change spectac- 
ularly between different nuclear spin states and can be dramat- 
ically affected by applied electric fields 10 . However, indirect 
methods are so far restricted to alkali-metal dimers and it will 
be challenging to extend them to other regions of the periodic 
table 15 . 

Direct cooling methods can be applied to a much larger va- 
riety of chemically interesting molecules, including OH, NH3, 
CO and LiH— ^. Stark deceleration, pioneered by Meijer 
and coworkers 16 , can be applied to polar molecules with large 
Stark effects, while helium buffer-gas cooling 20 has been par- 
ticularly successful for paramagnetic species. However, the 
temperatures so far achieved with direct methods are limited to 
tens of millikelvin, which is not cold enough to achieve quan- 
tum degeneracy. The development of a second-stage cooling 
method for such molecules is the biggest current challenge in 
the field. One of the most promising proposals is to use sym- 
pathetic cooling, which is based on the conceptually simple 
idea of bringing cold molecules into thermal contact with a 
bath containing ultracold atoms. So far sympathetic cooling 
has been successfully realized for ions 21,22 and some neutral 
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but not for neutral molecules. 



Linear molecules in spatially degenerate electronic states 
(n, A, etc.) are particularly attractive for Stark deceleration, 
as they exhibit first-order Stark effects at moderate electric 
fields (in contrast to molecules in £ states, which exhibit only 
second-order Stark effects). After deceleration, the molecules 
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can be loaded into traps where they are confined by static elec- 
tric or magnetic fields. Such static traps are not the only way 
to confine cold molecular species 25 , but they are experimen- 
tally the most accessible. In addition, atoms in open-shell S 
states (such as alkali-metal atoms, H( 2 S), N( 4 S), He( 3 S), and 
Cr( 7 S)), can be held in magnetic traps and may be suitable as 
coolants. 

Trapping with a static field is possible only if the atom or 
molecule is in a low-field-seeking state. However, the abso- 
lute ground state is always high-field-seeking. Thus, in addi- 
tion to the elastic collisions that lead to thermalization of the 
sample, there is always the possibility of inelastic collisions 
that transfer the colliding partners to a lower state and release 
kinetic energy. Inelastic collisions eject molecules from the 
trap and may lead to the heating of the sample. The success 
of sympathetic cooling therefore depends on the ratio of elas- 
tic to inelastic events, which should preferably be as large as 
possible. 

Molecular sympathetic cooling was first suggested for 
Rb+NH( 3 £~) 26 . Subsequently, potential energy surfaces 
and the appropriate collision cross sections have been 
calculated for a variety of candidate systems, includ- 
ing Mg+NH( 3 S-)^, Li+LiH( 1 £+)^22 Rb+NH 3 ^ and 
He+CH 2 ( 3 Si) 31 . Rb+ND 3 has also been explored experi- 
mentally 32 , though the inelastic collision rate in an electric 
field turned out to be too high for cooling. Studies of cold 
collisions with linear molecules in a II state in the presence of 
external fields have mostly been limited to cases when the sec- 
ond colliding partner is closed-shell. In particular, Tscherbul 
et al. 33 have investigated OH+He collisions and have shown 
how the inelastic cross sections can be reduced by combin- 
ing electric and magnetic fields to eliminate certain inelastic 
channels. Collisions of rotationally excited OH with He in the 
presence of electromagnetic field were analyzed by Pavlovic 
et al.—, while Bohn and coworkers 35 studied cold collisions 
between two OH molecules with long-range dipole-dipole in- 
teractions and concluded that the evaporative cooling of OH 
would be challenging. Lara et al— carried out theoretical 
studies of cold collisions of OH with Rb, taking account of 
multiple potential energy surfaces and including the hyperfine 
structure of OH. However, they did not include external field 
effects. 

There is thus a need for rigorous quantum studies of colli- 
sions between a II-state molecule and an open-shell S-state 
atom in the presence of external fields. In this paper, we 
extend the theory presented in Refs. 36 ' 37 to handle this case. 
This theory will be applicable to a broad set of experimen- 
tally important systems, including interactions of molecules 
such as OH, NO, CIO, and CH with alkali-metal and other 
magnetically trappable atoms. As an example, we present 
numerical results for collisions between OH( 2 n) and N( 4 S) 
in a magnetic field, with both colliding species initially in 



their fully spin-stretched low-field-seeking states. OH was 
one of the first molecules to be decelerated and trapped^ 3 ^, 
and many pioneering experiments with it have been reported. 
Gilijamse et al. 19 carried out a crossed-beam experiment, col- 
liding velocity-controlled OH molecules with Xe atoms; they 
were able to resolve state-to-state inelastic cross sections as a 
function of the collision energy. Similar experiments with im- 
proved sensitivity have recently been performed for OH col- 



liding with Ar, He, and D 
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An experiment to collide 



two velocity-controlled beams, of OH and NO, is in prepara- 
tion 42 . Sawyer et al. 38 have measured energy-dependent cross 
sections for collisions between magnetically trapped OH and 
slow D2 molecules. 

Tscherbul et al.— have recently suggested that spin- 
polarized nitrogen atoms are a promising coolant for sympa- 
thetic cooling experiments. N atoms at T >1 mK are stable 
against collisional relaxation between different Zeeman lev- 
els for a wide range of magnetic field strengths. Moreover, 
the low polarizability of the N atom leads to potential energy 
surfaces with an anisotropy much smaller than is usually en- 
countered for interactions with alkali-metal atoms. Theoret- 
ical and experimental studies for collisions of magnetically 
trapped N( 4 S) and NH( 3 £~) have been reported 4445 , showing 
that the trap loss in this system is fairly small and is caused 
mostly by the anisotropic magnetic dipole-dipole interaction 
between the atomic and molecular spins. 

This paper is organized as follows. In Sec. II we describe 
calculations of the high-spin (quintet) potential energy sur- 
faces resulting from interaction of the N( 4 S) atom with the 
OH( 2 n) molecule. In Sec. Ill we describe the effective Hamil- 
tonian used in the dynamical calculations and give the expres- 
sions for the matrix elements of the Hamiltonian. In Sec. IV 
we discuss the results of the scattering calculations and their 
implications for sympathetic cooling of OH by N atoms. Fi- 
nally, Sec. V summarizes and concludes the paper. 



2 Potential energy surfaces 

The interaction between the N( 4 S) atom and the OH( 2 n) 
molecule occurs on four adiabatic surfaces: 3 A', 3 A", 5 A', 
and 5 A". The triplet surfaces have been studied extensively 
to investigate the reaction N+OH— >NO+H that can take place 
on the 3 A" surface^ 5 ^. This reaction is the major source of 
the NO radical in the interstellar medium and is one of the 
key elementary processes in nitrogen chemistry. Formation of 
NO is barrierless, via a stable intermediate complex NOH, and 
is highly exothermic with 1.83 eV energy release. The other 
possible reaction channel N+OH— >NH+0 is energetically for- 
bidden for low-energy collisions. If we neglect minor spin- 
orbit coupling effects between the triplet and quintet states, 
the quintet surfaces are non-reactive. To our knowledge, the 
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Fig. 1 Contour plots of the quintet interaction potentials for N+OH: 5 A' (left-hand panel) and A" (right-hand panel). Energies are in cm 



quintet surfaces of N+OH have not been reported in the liter- 
ature thus far. 

We have carried out calculations of the quintet sur- 
faces using the unrestricted version of the coupled-cluster 
method with single, double, and noniterative triple excitations 
[UCCSD(T)]. The unrestricted version was chosen to circum- 
vent the problem of the lack of size-consistency for the in- 
teraction between two open-shell systems in spin-restricted 
coupled-cluster calculations^. The highly accurate aug-cc- 
pV5Z basis set of Dunning 53 was employed for all atoms and 
the counterpoise procedure^ was used to correct the com- 
puted interaction energies for basis-set superposition error. 
The MOLPRO suite of codes 55 was used in the electronic struc- 
ture calculations. 

Both the 5 A' and 5 A" potential energy surfaces were com- 
puted on a grid of points in Jacobi coordinates (R, 9), where 
R is the intermolecular distance measured from the centre of 
mass of 16 OH to the 14 N atom and 9 is the angle between the 
vector pointing from O to H in the OH molecule and the vec- 



Table 1 Characteristic points of the interaction potentials for the 
quintet states of N( 4 S) + OH( 2 IT). 



Rl ao 81 degrees V I cm 



Surface 



Global minimum 
Local minimum 
Saddle point 
Saddle point 



6.55 
6.36 
6.56 
6.66 



0.0° 
180.0° 

97.2° 
100.1° 



-71.5 
-61.0 

-45.8 



5 A', 5 A'' 



5 A' 
5 A" 



tor pointing from the centre of mass of OH to the N atom. The 
angle 9 — 0° thus corresponds to the linear O-H — N arrange- 
ment. The distance R was varied from 4.0 to 12.0 ao with an 
interval of 0.5 ao and from 12.0 to 20.0 ao with an interval of 
1 .0 ao . The angular grid points was chosen as the set of points 
for 11 -point Gauss-Lobatto quadrature, which include points 
at 9 = and 180°. The OH bond length was kept fixed at the 
monomer equilibrium value of 1.834 ao. 

Contour plots of the 5 A' and 5 A" potential energy surfaces 
are shown in Fig. 1. The shapes of the two quintet potentials 
are quite similar. The global minima appear for the linear ge- 
ometry O-H — N and have a depth of 120.9 cm -1 . There are 
also local minima 71.5 cm -1 deep, which occur at the linear 
N — O-H configuration. Note that for linear geometries the 
5 A' and 5 A" states are degenerate, so these minima are com- 
mon to the two surfaces. The set of stationary points of the 
potentials is completed by saddle points between the two min- 
ima, which are located at slightly different positions for the 
5 A' and 5 A" states. Table 1 gives the positions of the station- 
ary points on the surfaces and the corresponding interaction 
energies. The shapes of the quintet potential energy surfaces 
for N+OH closely resemble the high-spin (quartet) surface for 
N+NH reported by Zuchowski and Hutson 44 , although the 
global minimum for N+OH is about 30 cm -1 deeper than for 
N+NH. 

To perform quantum scattering calculations, it is necessary 
to expand the 5 A' and 5 A" surfaces in terms of angular func- 
tions. We adopt the convention of Alexander 37 and use spheri- 
cal harmonics in the Racah normalization Ck, q (9, <fi) (with an- 
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gle = 0) for angular representation of the potential. For the 
interaction of an S-state atom with a Il-state molecule, there 
are nonvanishing terms with q = and q = 2. The sum of the 
5 A' and 5 A" potentials is expanded in terms of functions with 
fc = 0, 



\[Va'(R,< 



v A »{R,e)\ = Y l c ki0 (o,o)v ko (R), (i) 



fc=0 



while the difference between the 5 A' and 5 A" potentials is 
expanded in terms of functions with k = 2, 



^[V A >(R,6)-V A „(R, 



= J2c k , 2 (S,0Wk2(R)- (2) 



fc=2 



Note that the definition of the difference potential, either 
Va' — Va" or Va" — Va>i depends in general on the sym- 
metry of the electronic wave functions of the interacting sub- 
systems 5657 . The radial functions Vk q {R) are obtained by 
projecting the sum or difference onto the appropriate angu- 
lar function, using Gauss-Lobatto quadrature to perform the 
numerical integration. Prior to this projection, the interpola- 
tion to obtain Va'(R,0) and Va"(R, 9) at an arbitrary value 
of R is done for each value of 6 using the reproducing ker- 
nel Hilbert space (RKHS) procedure 58 . For the quintet states 
of N( 4 S)+OH( 2 II), the dominant anisotropic term in the ex- 
pansion (Q]i is V2o(R) with a well depth of approximately 28 
cm -1 , while the dominant term in the expansion (0 comes 
fromV 22 {R). 

To improve the description of the potential at large R, we 
use an analytic representation in this region. Each radial com- 
ponent Vkq(R) is expanded at long range in terms of Van der 
Waals coefficients, 



V kq (R) = -^2C^/R n , 



(3) 



n— 6 



The expressions for the C^ q coefficients have been given by 
Skomorowski and Moszynski 57 , though with a different nor- 
malisation for q > from the one used here. We calculated 
the Van der Waals constants up to and including n=8, using 
the method described in Ref. 57 . The results are listed in Table 



Table 2 Long-range coefficients (in atomic units) for 
N( 4 S)+OH( 2 IT). 



(i 



k -t 

'6 
'6 







1 



27.84 



583.34 



51.60 



4.92 






1.23 


24.61 

-6.38 




312.00 




48.29 


159.09 




31.42 



2. For a weakly polarizable system such as N+OH, the neglect 
of higher-order coefficients with n > 8 is fully justified. We 
used the switching function of Janssen et al. 52 , with parame- 
ters a — 15 gq and 6 = 25 an, to join the asymptotic form 
based on the long-range coefficients and the RKHS interpola- 
tion of the ab initio points. 

3 Collision Hamiltonian 

3.1 Effective Hamiltonian 

We consider the case of an atom A( 2si+1 S), interacting with 
a diatomic molecule BC( 2 n), in the presence of an external 
magnetic field B. The direction of the field defines the labo- 
ratory (space-fixed) Z-axis. The system A-BC is described in 
Jacobi coordinates, with the r vector connecting the heavier 
and lighter of the atoms B and C, and R connecting the cen- 
tre of mass of BC and the atom A. By convention, lower-case 
and capital letters are used to represent the quantum numbers 
of the monomers and of the complex as a whole, respectively. 
The subscripts 1 and 2 refer to the monomers A and BC, re- 
spectively. For simplicity, the diatom will be treated as a rigid 
rotor in vibrational state v, although generalization to include 
its vibrations is straightforward. 

The Hamiltonian describing the nuclear motions of A+BC 
in the presence of magnetic field B can be written 



H = R 

2[i 



■i* R 

d& R 



L 2 



2nR 2 



Hr, 



Hi 



(4) 



where L is the space-fixed angular momentum operator de- 
scribing the end-over-end rotation of A and BC about one 
another and /i is the reduced mass of the complex. H mon 
contains all terms describing the isolated monomers, i.e. 
H mon = Hi + H-2- H\2 describes the interaction between 
the monomers: 



H\ 2 = H s 



V(R,t 



(5) 



Here, H SlS2 accounts for the direct dipolar interaction between 
the magnetic moments due to the unpaired electrons of the 
monomers, and V is the intermolecular interaction potential. 

If si 7^ and hyperfine terms are neglected, the Hamil- 
tonian for an isolated atom in the state 2si+1 S is fully deter- 
mined by the Zeeman interaction between the electron spin 
and the external magnetic field, 



Hi = gsUB Si • B, 



(6) 



where g$ is the electron g-factor, /ib is the electron Bohr mag- 
neton, and §i is the spin operator. 

The analogous Hamiltonian for a 2 n molecule can be writ- 



ten 



59 



H 2 — H ISO + Hz y2 + H\, 



(7) 
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where the rotational and spin-orbit contributions within the II 
state are collapsed into the first term, 



H r; 



B v n + A v I ■ s 2 . 



(8) 



B v and A v are the molecular rotational and spin-orbit con- 
stants, respectively, and n is the operator of the mechanical 
rotation of BC, which can be expressed as j — I — §2, where % 
I and §2 are the operators for the rotational, electronic orbital 
and spin angular momenta, respectively. H rso can be rewritten 

-ffrso = (A v +2B v )l z §2z 

+B V \f + P + $1 - 2j • s 2 + l 2 z - 2% l z ] . (9) 



The terms I , s\ and I 2 simply shift all the levels by a constant 
amount and are omitted below. The term H\, responsible for 
the A-doubling of the rotational levels of BC, is represented 
by the effective Hamiltonian 

Hx=J2 e ~ 2[q,Pr HvT 2 2q (j,j) + ( Pv + 2q v )T 2 2q (j,s 2 )} , 

q=±l 

(10) 
where <\> r is the azimuthal angle associated with the electron 
orbital angular momentum about the molecular axis defined 
by r, while p v and q v are empirical parameters. In Eq. dlOl . 
the second-rank tensor T 2 that couples two vectors ki and k2 
is defined as 



^(ki.ka) 



E< 

91,92 



l,gi;l,g 2 |2, g )Tj 1 (ki)Tj !1 (k 2 ), (11) 



where (l,gi; 1, 92 12, q) is a Clebsch-Gordan coefficient and 
the first-rank tensor components are Tj(k) = k z and 
T^_ 1 (k) = =F(kx i iky)/\/2. If only the electron spin and 
orbital contributions are taken into account, the Zeeman term 



Hz,2 = gs^B h ■ B + g' L [i B l ■ B, 



(12) 



where g' L is the orbital g-factor. For diatomic molecules of 
multiplicity higher than 2 (for example 3 II), an additional term 
describing the intramolecular spin-spin interaction must be in- 
cluded in the monomer Hamiltonian Q. 

The spin-spin dipolar interaction can conveniently be writ- 



ten 



59. 



H S1S2 = -9WbM^)V6 V(-l)«T2(s 1; s 2 )T 2 _ g (C) 



9 



(13) 

with T§((7) = C 2 , q (9, (j))FT z , where C 2 , q (9, <p) is a spherical 
harmonic function in the Racah normalization and (R, 6, (j)) 
is the set of relative spherical coordinates of the 'composite' 
atomic and diatomic electronic spins in the space-fixed frame. 
fio is the magnetic permeability of the vacuum. 



3.2 Basis sets and matrix elements 

The state of the BC molecule can conveniently be described 
using Hund's case (a) basis functions | A; S2&2; jurrij}, where 
S2 is the electron spin with projection 02 onto the molecular 
axis (body-fixed z axis), A is the (signed) projection of the 
electronic orbital angular momentum onto the molecular axis, 
and j is the angular momentum of BC with projections to onto 
the molecular axis and rrij onto the space-fixed Z-axis. For 
the body-fixed projections we have to = A + a. The state 
of the atom is characterized by the electronic spin function 
\sim Sl ). The basis set used here for the A-BC collision sys- 
tem is constructed as |sito Si ) |A; S2CT2; jurnij) \LMl), where 
\LMl) are functions describing the relative motion of A and 
BC in the space-fixed reference frame. 

In the presence of a magnetic field, the conserved quanti- 
ties are the projection M tot of the total angular momentum, 
Mtot = rris! + rrij + Ml, and the total parity V. An electric 
field would mix states of different total parity. In the absence 
of an electric field it is most efficient to use a parity-adapted 
basis set, \sim Sl ) \s2', juirrije) \LMl), with 



|s 2 ;jwm i e) 



|l;s 2 <7 2 ; jujrrij 



1 

; 71 



+ e{-iy- S2 \-l; S2 -V2;j-um 



31 ' 



(14) 



where u) = \u\, 02 = U) — 1 and e = ±1. In this basis set, 
the parity of BC is P2 = e(—iy~ S2 , and that of the triatomic 
system is V = p\P2{— 1) L - The matrix elements of L 2 and 
H\ are diagonal, and given by h 2 L(L + 1) and gs^m Sl B, 
respectively. 

We next give the matrix elements of all terms in the Hamil- 
tonian of Eq. ©, although only those involving the atomic 
spin are new in the present work. The terms that do not in- 
volve atomic spin are the same as for collisions with a closed- 
shell atom and were previously given by Tscherbul et a/. 33 . 
However, the published version of Ref. 33 contains a number 
of typographical errors, so we report the correct expressions 
here. 

The matrix elements of the molecular rotation/spin-orbit 
operator are 

(s2;juimje\ H rso |s 2 ;ja>'mje) 

<W{(4, + 2B v )(6j-1) + B v [j(j + 1) - 2u 2 } } 

-B nl <W-ia-0',w')a_(s2,w' - 1) 

+8uu> + ia + (j,u}')a + (s2,u}' - 1) , (15) 



where we use a±(j, m) = \/j{j + 1) — m(m ± 1) both to 
simplify the equations and to ease comparison with Ref. 33 . 
The off-diagonal terms on the right-hand side connect differ- 
ent spin-orbit manifolds related by uj' = u> ± 1. 
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The A-doubling matrix elements are 

{s 2 ;ju)mje\H\ \s 2 ; jui'mjc) = -e(-l) J ~ B2 a-(j,ui') 

x [5Q2-ui>q v a-(j,u)' - 1) - foi-o'O^ + %Qv)oi+(s2,u' - 1)] 

(16) 

and also couple states with different uj. For a 2 IT molecule, the 
first factor inside the square brackets mixes the 1/2 and 3/2 
states, while the second is non-zero only for u = ui' = 1/2. 

The matrix elements of the Zeeman interaction for the 
molecule BC are 



(s 2 \ jumj(\ Hia \s 2 ;j'uj'm j e) 

= MB B(-ir^'[(2j + l)(2/ + l)] 1/2 

a+(s 2 ,ui' - 1) ( j 1 



ffs- 

3 

Co 



V2 



gs- 



-(82,Q'-1) 



V2 



+ \g s (ui-l)+g' L \ 



3 1 

uj 



(17) 



and mix both different rotational and different spin-orbit 
states. 

To determine the matrix elements of the spin-spin dipo- 
lar interaction it is natural to expand the second-rank ten- 
sor T 2 (Ji,S2) as a linear combination of the products of 
the space-fixed components of first-rank tensors T* (si) and 
Ti 2 (s2). The matrix elements of T^fsi) can be calculated 
directly in our basis set, while for Tp 2 (,§2) we first need to 
transform from the space- to the body-fixed frame, 



Tp 2 (s2)=E^( n ) T 3^). 



(18) 



where T> J KM is a Wigner rotation matrix and il represents the 
Euler angles for the transformation. The matrix elements in 
the primitive basis set are 

(LM L ; A; s 2 a 2 ; jurnif, Sim sl \H BlB2 \sim' ai ; A; s 2 a' 2 ;j'uj'm'j- V ' M' L ) 
^-VS0\ slS2 {R)(-l) sl - m ^ +S2 ~' 72+m ^ u '~ ML 
x [si(si + l)s 2 (s 2 + l)(2si + l)(2s 2 + l)(2j + l)(2j' + 1)] 1/2 

,i/2 f L 2 U 

I ^000 

1 1 2 W si 1 si 

-P ) V -m ai pi m' Sl 



x [(2L + 1)(2L' + 1) 

x E 



Pi P2 
1,P2,9 " 

S 2 1 S2 

-<T2 q v 2 

3 1 / 

-lo q uj 



3 1 / 

— rrij p 2 m'j 

L 2 l! 

-M L -p M' l 



and the corresponding matrix elements in the parity-adapted 
basis set are 



{LM L ;s 2 ;jum j e; Sim 31 1 H Sl s 2 \sim' si ; s 2 ; j'uj'm'jC; L' M L ) 

= Vs6\ sl s 2 (R)(-lY 1 ~ msi+S2+m3+2 ' : " ML 

x [si(si + l)(2si + l)s 2 (s 2 + l)(2s 2 + l)(2j + l)(2j' + 1)] 1/2 



[(2L + 1)(2L' + 1) 



1/2 



L 2 L' 




1 1 

Pi P2 



-P 



- E 

P1.P2.9 

S 2 1 S2 

-cj + 1 q uj' — 1 

i i nf £ 

-uj g uj 1 I \ —Ml 



Si 1 si 

-m 31 pi m' si 



x 



X 



.7 



1 3 
—rrij p 2 m'j 

2 L' 
-P M' L 



(20) 



where p = p\ + p 2 , A S1S2 (-R) = £^aoa 2 /-R 3 is the R- 
dependent spin-spin dipolar coupling constant and a s=y 1/137 
is the fine-structure constant. 

Finally, the matrix elements of the interaction potential are 



(LM L -S2;jojmje\V\s 2 -,j'uj'm'je'-,L'M' L ) 



(-1)' 



,-— Q 1 — m 



M(2j + 1)(2/ + 1)(2L + 1)(2L' + 1)] 1/2 



Y^\[l + ee'{-l) k ]{-l) 



k.nik 



'"A, [ ■? fe 

TOj 771 fc 



L k L' 



3 k f 

uj -uj' 



L k V 

-M L m k M' L 



V k0 (R) 



(1-<WK 



3 k j' 

uj -2 uj' 



V k2 (R) 
(21) 



(19) 



where V k o (R) and V k2 (R) are the radial strength functions of 
Eqs. (fl} and (f2]). It is readily seen that states belonging to the 
same spin-orbit manifold are coupled through the 'average' 
of the A' and A" potential surfaces, while those of different 
manifolds are connected through their difference. In addition, 
the factor 3 [l + ee'(— l) k ] guarantees that states of the same 
monomer parity are connected by terms V kq (R) with even k, 
while those with odd k couple rotational levels of opposite 
parity. It follows from this that a strong parity-conserving 
propensity rule for transitions involving different spin-orbit 
manifolds can be expected. 
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4 Dynamical calculations 

4.1 Computational details 

Expanding the Schrodinger equation with the Hamiltonian of 
Eq. (|4|i in the parity-adapted basis set ( TT4T > yields a set of 
coupled differential equations. We have written a plug-in for 
the MOLSCAT general-purpose quantum molecular scatter- 
ing package— to implement the matrix elements described 
above for collisions between an open-shell S-state atom and a 
2 II-state molecule in a magnetic field. We solved the coupled 
equations numerically using the hybrid propagator of Alexan- 
der and Manolopoulos 61 , propagating from i? m i n = 4 a to 
Rmid = 25 a using a fixed-step log-derivative propagator 
with interval size 0.02 a,Q and from i? m id to i? max = 800 a$ 
using a variable-step log-derivative propagator based on Airy 
functions. MOLSCAT applies scattering boundary condi- 
tions at i? m ax to extract scattering S-matrices, which are then 
used to calculate elastic and inelastic cross sections. 

Values of the OH molecular constants in the monomer 
Hamiltonian were taken from Refs. 62f)3 . After performing nu- 
merous test calculations, we decided to include basis functions 
withj < 9/2 and L < 8, which gives convergence of the cross 
sections to within approximately 1%. 

4.2 Results 

The lowest rotational state of OH in its ground X 2 H state at 
zero field is a A doublet with j = 3/2. The doublet con- 
sists of two states, referred to as e and /, which have opposite 
parity and are separated by 0.059 cm -1 , with \j = 3/2, e) 
being the ground state. A magnetic field splits each compo- 
nent of the doublet into four states differing by the projec- 
tion of the angular momentum riij on the field axis (rrij = 
3/2, 1/2, -1/2, -3/2). For the N atom in its 4 S ground state, 
a magnetic field produces four Zeeman levels, with spin pro- 
jections m Sl = 3/2, 1/2, —1/2 and —3/2. The combination 
of 8 Zeeman levels of OH with 4 of the N atom yields 32 
asymptotic levels (thresholds), as shown in Fig. [2] In princi- 
ple, even at zero field each of the levels is further split due to 
hyperfine interactions, although in the present work hyperfine 
effects are neglected for simplicity. 

N and OH can both be magnetically trapped in their spin- 
stretched states, with m si = 3/2 and rrij — 3/2, respec- 
tively. There are two such states for OH, originating from the 
e and / components of the A doublet. We choose the initial 
state to be \m Sl = 3/2)|mj = 3/2, e), shown with a red 
line in Fig. |2] This is likely to be more favourable for sympa- 
thetic cooling than \m Sl — 3/2)\rrij — 3/2, f) (shown with 
a solid blue line in Fig. |2), because there are fewer inelastic 
channels open for Zeeman relaxation at low collision ener- 
gies. In particular, we avoid transitions between the two fully 
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Fig. 2 Energy levels of noninteracting N( 4 S)+0H(X 2 I1, j = 3/2) 
in a magnetic field. The solid red and blue lines indicates the 
spin-stretched low-field-seeking states \m si = 3/2}|rrij = 3/2, e) 
(red) and \m 31 = 3/2) \mj = 3/2, /) (blue). The dotted blue line 
shows state \m Sl — 3/2)\nij — 1/2, /). 



\m Sl = 3/2) \rrij = 3/2, e), at collision energies below about 
85 mK. The only centrifugal suppression in such a process, 
even for an incoming s wave (Li = 0, Af^.i = 0) is due to a 
p-wave barrier in the outgoing channel (Lf — 1, Mlj — 0) 
with a height of only 1 1 mK, necessitated by the change in OH 
monomer parity. 

The interaction between collision partners that are initially 
in fully spin-stretched states takes place almost entirely on 
the quintet (high-spin) potential energy surfaces. A full de- 
scription of exit channels in which m Sl + rrij has changed re- 
quires triplet surfaces, but including these explicitly would be 
computationally prohibitively expensive. In the present work, 
we effectively approximate the triplet potential surfaces with 
the corresponding quintet ones. This approximation is closely 
analogous to that used for N+NH in ref. 44 . 

In a low-energy inelastic collision, the quantum state of 
at least one of the colliding species changes and kinetic en- 
ergy is released. There are two main mechanisms that pro- 
duce inelasticity in ultracold collisions of an open-shell S- 
state atom with a molecule in a 2 II state. The first is direct 
coupling through the anisotropy of the interaction potential, 
which drives transitions to states with the molecular quantum 
number nij reduced by at least 1 and the atomic spin projec- 



spin-stretched states, from \m Sl = 3/2)|toj = 3/2, f) to tion m Sl unchanged. This mechanism is also present in col- 
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Fig. 3 Elastic and total inelastic cross sections for N+OH scattering 
at different magnetic field strengths B. The elastic cross section is 
almost unaffected by the field strength and is shown only for 
B = 10G. 



tion is almost unaffected by the magnetic field. At ultralow 
collision energies, the inelastic cross sections are suppressed 
due to centrifugal barriers in the outgoing channels, and the 
total inelastic cross section grows with increasing field be- 
cause the increasing kinetic energy release helps overcome 
these barriers. For example, the energy released by relaxation 
to the state \m Sl = 3/2)|toj = -1/2, e) at a field of 560 G 
is sufficient to overcome the cf-wave barrier, and Fig.[5]shows 
how the inelastic cross section is enhanced for fields of 500 
G and higher in the s-wave regime. For collision energies be- 
tween 4 mK and 80 mK, both the elastic and inelastic cross 
sections feature numerous resonances, mostly Feshbach reso- 
nances due to coupling with higher-energy closed channels. 

As discussed above, there are two mechanisms driving 
transitions between different Zeeman levels, one driven by 
the spin-spin dipolar term -ffsis 2 ar, d the other driven by the 
anisotropy of the intermolecular potential Vkq (R). The mech- 
anism involving H Sl s 2 dominates for low fields (10 G and be- 
low) and in the s-wave regime. For higher fields (100 G and 



lisions between a closed-shell atom and a 2 II molecule and 
has been described by Tscherbul et al.—. The second mech- 
anism arises from coupling by the spin-spin dipolar interac- 
tion H Sl s 2 - Here, the final Zeeman state may have quantum 
numbers rrij and m Sl reduced by at most one. Such pro- 
cesses are also present in collisions of an open-shell S-state 
atom with 2 S or 3 S molecules, or indeed between two alkali- 
metal atoms. Collisions of spin-polarized S-state atoms with 
2 II molecules thus combine two direct mechanisms for cou- 
pling between different Zeeman levels. 

The most important contribution to coupling by the interac- 
tion potential comes from the anisotropic term V2o(R), which 
induces direct transitions from the OH state \irij — 3/2, e) 
to \rrij — 1/2, e) and \irij = — 1/2, e). This occurs even in 
the s-wave regime (Li = 0). An s-wave collision in which 
rn Sl + rrij decreases requires Mlj > 1 to conserve Af to t. If 
the monomer parity is unchanged, conservation of total par- 
ity then requires Lf > 2. There is thus a centrifugal barrier in 
the outgoing channel, which suppresses the inelastic cross sec- 
tions for low collision energies and low fields. For N+OH, the 
centrifugal barriers are relatively high due to the low reduced 
mass and small Cg° coefficient: the height of the ri-wave bar- 
rier is 7 1 mK. 

Fig.[3]shows the cross sections for Zeeman relaxation in col- 
lisions of OH(X 2 n, \m,j = 3/2, e» with N( 4 S, \m si = 3/2)) 
for magnetic field strengths B = 10, 100, 500 and 1000 G. 
At low collision energies (below 0. 1 mK), the cross sections 
behave according to the Wigner threshold laws 64 : the elastic 
cross section is constant, while the total inelastic cross section 
grows with decreasing energy as E^ 1 / 2 . The elastic cross sec- 
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Fig. 4 Comparison of the total inelastic cross sections (upper panel) 
and the s-wave contributions to them (lower panel) for N+OH, 
obtained with the spin-spin dipolar interaction included or neglected 
in the Hamiltonian, for magnetic fields B = 10 and 100 G. 
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Fig. 5 Comparison of the s wave total inelastic cross sections for N+OH with those obtained with either the spin-spin dipolar term or the 
anisotropy of the interaction potential neglected. Left-hand panel: B — 10 G; right-hand panel: B — 100 G. 



above), the opposite is true and the relaxation is driven by 
Vkq(R). Fig. [4] shows the integral cross sections and the s- 
wave contribution for the two lowest fields (10 and 100 G), 
with the H SlS2 term in the Hamiltonian included or neglected. 
Fig. [5] compares the s-wave contributions for the same two 
fields with those obtained by neglecting either the spin-spin 
dipolar term -ff SlS2 or all the anisotropic terms Vk q (R). At 10 
G, -ff S is 2 greatly enhances inelastic processes in the ultracold 
regime: at 10 -5 K, the enhancement is almost two orders of 
magnitude. However, for B = 100 G, Vk q (R) is dominant 
over the whole range of energies. 

The way in which the spin-spin dipolar interaction induces 
transitions between different Zeeman levels is exactly par- 
allel to that described by Janssen et a/. 65 '- 6 . It is a purely 
long-range effect caused by narrowly avoided crossings be- 
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N( 4 S) OH( 2 rU = 3/2,e) 
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Fig. 6 Pattern of first-order couplings between different Zeeman 
levels of N( 4 S)+OH(X 2 II, j = 3/2) through the spin-spin dipolar 
interaction and the anisotropy of the interaction potential, for 
incoming s wave (Li — 0) and outgoing d wave (L/ = 2). 



tween the potential adiabats at very long range, which enable 
transitions between Zeeman levels without the need to pene- 
trate centrifugal barriers. In the present case, avoided cross- 
ings due to the dipolar term are present between the adia- 
bat asymptotically correlating with the incident s-wave chan- 
nel \m Sl = 3/2) \rrij = 3/2, e) and other adiabats corre- 
lating with the states \m Sl = 3/2) |mj = 1/2, e), \m Sl — 
1/2) |mj = 3/2, e), and \m Sl = l/2)|mj = 1/2, e). The 
p-wave and higher-L contributions to the total inelastic cross 
sections are almost unaffected by the inclusion of -ff S is 2 f° r 
any field and collision energy. This arises because the long- 
range avoided crossings for incident channels with centrifugal 
barriers are energetically inaccessible at low energies. 

Channels corresponding to different Zeeman levels are also 
directly coupled by the anisotropy of the intermolecular po- 
tential Vkq(R). Fig. [6] shows a schematic illustration of the 
first-order couplings by -ff S is 2 an d Vk q (R) for collisions in- 
volving an incoming s wave and outgoing d waves. Because 
of this, long-range avoided crossings are present even if we 
neglect H SlS2 . However, the effect of the avoided crossings 
on collision outcomes is much more pronounced for crossings 
due to H Sl s 2 man f° r those due to Vk q (R). The latter dies off 
much faster with R than H SlS2 , and is one or two orders of 
magnitude weaker at the positions of the long-range avoided 
crossings. The ratio of the coupling strengths is approximately 
H SlS2 (R)/V 2 o(R) - E h a 2 /Cl°(R/a )- 3 . The avoided 
crossings for B — 10 G occur at distances ranging from 159 
to 342 an, corresponding to a ratio -ff S is 2 (R)/V2o(R) between 
10 and 100. It follows from an approximate Landau-Zener 
model 67 that the probability of ending in a different asymp- 
totic level after a nonadiabatic transition is proportional to the 
square of the coupling between the diabats if the coupling is 
relatively small. 
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Fig. 7 State-to-state inelastic cross sections (s-wave contribution 
only) for fields of 10 G (upper panel) and 100 G (lower panel). 



The interplay between the spin-spin dipolar term and the 
intermolecular potential terms is also manifested in the state- 
to-state cross sections. Fig. Q shows state-to-state cross sec- 
tions (s-wave contributions only) for B — 10 G and 100 G. 
At B = 10 G, in the region where the -ff S is 2 term dominates 
(below 1 mK), the most important transitions are to states with 
rrij or m Sl quantum numbers reduced by 1, which are those 
coupled to the incident channel \m Sl = 3/2)\rrij = 3/2, e) 
by H SlS2 , while for collision energies above 1 mK the dom- 
inant inelastic channels become |rn Sl = 3/2) \rrij — 1/2, e) 
and \m Sl = 3/2}\m,j = —1/2, e), which are those coupled by 
Vk q {R). At B = 100 G, only channels coupled by Vk q (R) are 
important. 

The s-wave cross sections at B = 10 G exhibit two distinct 
resonant structures: a strong feature near 15 mK and a weaker 
one around 41 mK. Both are Feshbach resonances caused by 
coupling to closed channels arising from the / component of 
the A doublet of OH. The coupling arises almost exclusively 
from the Vio {R) term in the intermolecular potential, which 
couples states of different monomer parity. The Feshbach res- 




Fig. 8 The lowest adiabatic potential energy curves for A/ to t = 3 
and I/max = 2 correlating with thresholds with the state of the N 
atom unchanged (\m si = 3/2)), at B — 10 G. Two solid horizontal 
lines indicates the position of the bound states responsible for the 
two sharp Feshbach resonances in the s and ri-wave contribution to 
the inelastic cross sections. 



onance near 15 mK can be attributed to a bound state on the p- 
wave adiabat correlating with the \m si = 3/2} jm^ = 1/2, /) 
threshold, as shown in Fig. [8] This resonance moves to smaller 
energies with increasing field, because the energy difference 
between the \rrij = 3/2, e) and \rrij = 1/2, /) states (red 
and dotted blue lines in Fig. [2] respectively) decreases as the 
field increases. For sufficiently large field (B > 1200 G), 
this resonance will disappear as the \m,j = 1/2, /) level drops 
below \m,j = 3/2, e). The second Feshbach resonance near 
41 mK can be attributed to a bound state on the p-wave adi- 
abat correlating with the \m Sl = 3/2) \rrij = 3/2, f) thresh- 
old. The position of this resonance is almost unaffected by 
the field strength since the energy difference between the two 
spin-stretched states, \rrij = 3/2, e) and \rrij — 3/2, /), is 
independent of magnetic field. 

The fact that these are Feshbach (rather than shape) res- 
onances is confirmed by several observations. First, the ri- 
wave contributions to the inelastic cross sections show reso- 
nant structures at exactly the same energies as the s-wave con- 
tribution. Secondly, the positions and shapes of the Feshbach 
resonances can be reproduced using even the smallest possi- 
ble basis set that allows inelastic transitions, with j < 3/2, 
L < 2, and potential terms Vk q (R), k < 2. Thirdly, the pres- 
ence of the T4o (R) term, which does not couple the incident 
and outgoing channels directly, is crucial for the existence of 
the resonances. It is worth noting that no such structure due to 
Feshbach resonances would be present for collisions involving 
the initial state \m Sl — 3/2)|mj = 3/2, f), with OH in the 
upper component of its A doublet, since no low-lying closed 
channels are present in that case. However, molecules in the 
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Fig. 9 Ratio of elastic to total inelastic cross sections for N+OH at 
different magnetic fields. 



/ state are likely to undergo fast relaxation to the e state in 
collisions driven directly by Vio(R). 

Fig. [9] shows the ratio of the elastic to total inelastic cross 
sections as a function of collision energy. The ratio is not 
favourable for sympathetic cooling of OH by collision with 
ultracold N atoms, except at fields below 10 G and collision 
energies below 1 mK. The cross sections presented here may 
be compared to those for N( 4 S)+NH( 3 £~) by Zuchowski and 
Hutson 44 . The ratio of the elastic to inelastic cross sections 
is at least an order of magnitude lower for N+OH than for 
N+NH. Two main reasons for this can be identified. First, 
the spin-stretched component of the rotational ground state of 
NH ( 3 £~, n = 0) is not directly coupled by the potential 
anisotropy to any other accessible Zeeman level, whereas such 
a coupling does exist for the ground state of OH( 2 n, j = 3/2) 
(or any other molecule with j > 1). Secondly, there are low- 
lying states arising from the / component of the A doublet 
in the OH radical that create many Feshbach resonances and 
increase the inelasticity. Both effects are particularly strong 
for collision energies above 10 mK, where the contributions 
from p and d incoming waves to the inelastic cross sections 
are dominant. For all field strengths, the ratio of elastic to 
inelastic cross sections at collision energies above 1 mK is 
more than 10 times larger for N+NH than for N+OH. 



4.3 Potential dependence 

The results of scattering calculations at ultralow collision en- 
ergies are in general very sensitive to the details of interaction 
potentials. To estimate the accuracy of the calculated poten- 
tial energy surfaces for N+OH, we have carried out additional 
electronic structure calculations for the geometry correspond- 
ing to the global minimum of the potentials at the linear N-OH 



geometry. In the aug-cc-pV5Z basis set, the global minimum 
has a well depth of 120.9 cm -1 , while in the aug-cc-pV6Z ba- 
sis set this shifts to 121.8 cm -1 . Based on these two results, 
we can estimate the complete basis set limit of the CCSD(T) 
method to be 122.9 cm -1 , using the extrapolation formula for 
correlation energy as given in Ref. 68 . This corresponds to an 
error estimate of 1.7% for our full potential surfaces using the 
aug-cc-pV5Z basis set. To estimate the error in the correla- 
tion energy obtained from the CCSD(T) method, we have per- 
formed full configuration-interaction (FCI) calculations with 
eight electrons correlated in the cc-pVDZ basis set. The rela- 
tive contribution of the FCI correction to the CCSD(T) result 
should only be weakly dependent on the basis set used, so even 
in the small cc-pVDZ basis set we should obtain a reliable es- 
timate of the FCI valence-valence correlation correction. The 
FCI correction to the CCSD(T) result accounts for approxi- 
mately 1.5% of the interaction energy at the global minimum. 
We can thus estimate the uncertainty of our potential energy 
surfaces to be 4% at worst. 
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Fig. 10 Cross sections obtained with the interaction potential scaled 
by a constant factor, V(R) — ¥ A ■ V(R), for magnetic field B — 10 
G and collision energies of 10 /iK (upper panel) and 1 mK (lower 
panel). 
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To assess the sensitivity of the scattering results to the un- 
certainty in the interaction potential, we have carried out cal- 
culations with the interaction potential scaled by a constant 
factor A in the range 0.96 < A < 1.04, corresponding to the 
estimated error bounds in the calculated potential energy sur- 
faces. The results at B = 10 G are shown in Fig. \W\ for 
collision energies of 1 mK and 10 /jK The weak dependence 
of the cross sections on the potential scaling is disturbed by 
the presence of sharp resonances, which occur when bound 
states of the N-OH complex cross the incoming threshold (or 
more precisely the collision energy) as a function of A. Two of 
the peaks in the inelastic cross sections, near A = 1.010 and 
A = 1.026, can be attributed to the Feshbach resonances in 
the s and d partial-wave contributions discussed above. The 
two additional peaks at A = 0.97 and A = 1.017, which 
broaden substantially with collision energy, are due to shape 
resonances in the p-wave partial cross section. If the true po- 
tential happens to bring one of these resonances close to zero 
energy, it may change the ratio of elastic to inelastic cross sec- 
tions quite dramatically. However, Fig.[l0]shows that the reso- 
nances occur in quite narrow ranges of A, so that there is a low 
probability that the true potential will be such that the ratio of 
elastic to inelastic cross sections is seriously affected by reso- 
nances for collision energies below 1 mK. It may also be noted 
that the numerical results obtained with the unsealed potential 
(A = 1) are fairly typical of the range expected for N+OH 
on plausible interaction potentials, in the sense that the low- 
energy elastic cross section (around 1000 A 2 ) is close to the 
value a = Ana 2 = 712 A 2 obtained from the mean scattering 
length a defined by Gribakin and Flambaum 69 . 

5 Summary and conclusions 

We have presented a theoretical study of the relaxation pro- 
cesses in collisions between an atom in an open-shell S state 
and a molecule in a 2 IT state, in a magnetic field, using the 
example of N( 4 S)+OH( 2 II). The transitions between dif- 
ferent Zeeman levels in such collisions are driven by two 
mechanisms: coupling through the spin-spin dipolar term 
and through the anisotropy of the interaction potential. Both 
mechanisms are present in first order. The spin-spin dipolar 
term dominates when both the collision energy and the mag- 
netic field are low, while the anisotropy of the interaction po- 
tential dominates at higher energies or fields. In the latter 
regime, the spin-spin dipolar term can be neglected. Neglect- 
ing the dipolar interaction is equivalent to treating the atom as 
closed-shell, which dramatically reduces the cost of the scat- 
tering calculations. 

An important general point is that spin relaxation collisions 
can be driven directly by the anisotropy of the interaction po- 
tential for any molecule that has rotational angular momen- 
tum. Since the anisotropics of atom-molecule and molecule- 



molecule interaction potentials are typically quite large, this 
will often provide an important trap loss mechanism for such 
states. For molecules in 2 II states, this is true even for the 
molecular ground state. 

For the case of N+OH, the spin-spin dipolar term dominates 
at collision energies below about 1 mK and magnetic fields of 
10 G or less. In this regime, the ratio of elastic to inelastic 
cross sections is greater than 100 and thus favourable for sym- 
pathetic cooling. However, if either the collision energy or the 
magnetic field is significantly above this, inelastic processes 
due to the potential anisotropy dominate and the ratio of elas- 
tic to inelastic cross sections falls. This suggests that sympa- 
thetic cooling of OH by collisions with N atoms is unlikely to 
be successful except at collision energies below 1 mK. 
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